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We investigate the orbital origin of the Fano-Kondo line shapes measured in STM spectroscopy 
of magnetic adatoms on metal substrates. To this end we calculate the low-bias tunnel spectra of 
a Co adatom on the (001) and (111) Cu surfaces with our density functional theory-based ab initio 
transport scheme augmented by local correlations. In order to associate different d-orbitals with 
different Fano line shapes we only correlate individual 3d-orbitals instead of the full Co Sd-shell. We 
find that Kondo peaks arising in different d-levels indeed give rise to different Fano features in the 
conductance spectra. Hence the shape of measured Fano features allows us to draw some conclusions 
about the orbital responsible for the Kondo resonance, although the actual shape is also influenced 
by temperature, effective interaction and charge fluctuations. Comparison with a simplified model 
shows that line shapes are mostly the result of interference between tunneling paths through the 
correlated d-orbital and the sp-type orbitals on the Co atom. Very importantly, the amplitudes of 
the Fano features vary strongly among orbitals, with the Sa^-orbital featuring by far the largest 
amplitude due to its strong direct coupling to the s-type conduction electrons. 


I. INTRODUCTION 

The Kondo effect is one of the most fascinating phe¬ 
nomena in condensed matter physics, occurring in a vast 
number of different systems (see, e.g.. Ref. [T] and ref¬ 
erences therein), racing from bulk metals doped with 
magnetic impurities^lHl to nano scale systems such as 
semiconductor quantum dotm and carbon nanotubes 
connected to metal lead^^. Generally, the Kondo effect 
leads to the quenching of a local magnetic moment asso¬ 
ciated with localized and strongly interacting electronic 
states in the system by interaction with the conduction 
electrons. The quenching of the spin is accompanied by 
drastic changes in the electronic and transport proper¬ 
ties. This strong impact on the electronic and magnetic 
properties of a system makes the Kondo effect an impor¬ 
tant factor for the functionality of atomic and molecular- 
scale electronic devices. 

Since the pioneering works of Li et and Madhavan 
et scanning tunneling spectroscopy (STS) has be¬ 
come a standard tool for probing the Kondo effect of mag- 
netic adatoms and molecules on metallic substrates^JS. 
The Kondo effect arises from the interaction of the mag¬ 
netic moment of the adsorbate with the conduction elec¬ 
trons of the metal surface, and leads to the screening of 
the magnetic moment by formation of a total spin singlet 
state with the conduction electrons. The formation of 
the Kondo-singlet state is signaled by the appearance of 
a strongly renormalized quasi-particle peak at the Fermi 
level, the so-called Abrikosov-Suhl or Kondo resonance. 
In STS the appearance of the Kondo peak in the local 
DOS of the atom or molecule leads to a zero-bias anomaly 
(ZBA) in the tunnel spectra which is generally well de¬ 
scribed by a Fano line shap^j^, although it has recently 
been found that the ZBAs are actually much better de¬ 
scribed in terms of generalized Frota line shape^^ as the 
Frota function yields a much better descri ption of the 
Kondo peak than the Lorentz function i^^l^^ l 

The origin of the Fano-like line shape is either un¬ 


derstood as due to the interference of different tunnel¬ 
ing paths - one via the strongly interacting orbitals of 
the magnetic atom bearing the sharp Kondo resonance, 
and others going directly to the substrat^^MUl _ or is ex¬ 
plained in terms of tunneling into the surface alon^SHESl 
A recent studjE^ combines density functional theory 
(DFT) with numerical renormalization group (NRG) cal¬ 
culations and determines the line shape by looking at 
energy-dependent transmission eigenvalues. Surprisingly, 
no systematic study of the relation between orbital sym¬ 
metry of the orbital(s) bearing the Kondo resonance and 
the shape of the resulting Fano resonances has been con¬ 
ducted so far. 

In this paper, we intend to close this gap by calculating 
the Fano line shapes corresponding to Kondo peaks ap¬ 
pearing in different orbitals of the Sd-shell of a magnetic 
atom on metal surfaces. To this end we select individ¬ 
ual d-orbitals and perform ab initio quantum transport 
calculations augmented by local correlations for the se¬ 
lected d-orbital only. There is merit in doing so: Even 
in a multi-orbital situation the Kondo effect is signaled 
by Kondo peaks in individual d-orbitals, and often the 
Fano-Kondo feature of one d-orbital will be dominant in 
the tunnel spectrum due to different tunneling matrix 
elements and Kondo scales. 

We choose to study Co@Gu(001) and Co@Cu(lll) as 
our test systems, which have b een extens ively studied 
theoretically and experimentalljIi^IIIIlIIlMlS, We find 
that Kondo peaks arising in different d-levels indeed give 
rise to different Fano features in the conductance spec¬ 
tra. However, temperature, effective interaction and oc¬ 
cupancy of the d-orbital also play an important role. 
With one notable exception, a simplified two-level model 
consisting of the d-orbital bearing the Kondo resonance 
and one s- or p-orbital on the adatom accounts for the 
calculated line shapes. This shows that in these cases 
tunneling into substrate states only plays a minor role 
for determining the actual line shapes. 

The paper is organized as follows: In Sec.|Tl]we briefly 
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core electrons (3s3p) of the Co and Cu atom^^lHlSl xhe 
effect of the bulk electrodes S and T, which are modeled 
by Bethe lattice^^, on the electronic structure of D is 
taken into account via self-energies Ss and Et- The KS 
Green’s function (GF) of the D region is thus given by: 
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FIG. 1. Schematic drawing of an STM tip probing a mag¬ 
netic atom on a metal substrate. The system is separated 
into three parts: the device region D (gray/yellow) contains 
the magnetic atom (dark gray/red) hosting the d-orbital giv¬ 
ing rise to the Kondo peak, and parts of the substrate and 
STM tip. T and S (light gray/blue) are the bulk electrodes 
connected to the STM tip and substrate, respectively. 


describe the method for calculating the zero-bias anoma¬ 
lies in the conductance spectra corresponding to Kondo 
peaks in different orbitals. In Sec. [m] we introduce two 
types of Fano line shapes: the standard one based on 
a Lorentzian resonance for the localized state and one 
based on the Frota line shape better suited for describing 
the Kondo resonance. In Sec. |IV| we present results for a 
Go adatom placed on a Gu(OOI) and a Gu(lll) surface, 
respectively. In Sec. we devise a simplified model cap¬ 
turing the essence of the different situations encountered 
for different orbital symmetries and discuss the obtained 
results in the context of this model. In Sec. IVIl a more 


general discussion follows relating our results to other 
experimental and theoretical works. Finally, in Sec. |VII[ 
we conclude this work with some general remarks on the 
significance of our results for other atomic or molecular 
Kondo systems. 


II. METHOD 


G°(cc) = ((a; + Ai)PD-i^S-STH-EsM) \ (1) 

where /i is the chemical potential, Pd the projection op¬ 
erator onto D and is the KS Hamiltonian of the D 
region. 

In order to capture Kondo physics, electronic correla¬ 
tion beyond conventional DFT have to be included. This 
is done by combining DFT with the one-crossing approxi¬ 
mation (O GAI^ , following the scheme developed in pre¬ 
vious worlPSl^. In contrast to previous work, we are 
interested in the Kondo signatures of specific d-orbitals, 
and not of the entire 3d-shell. 

Hence we add a Hubbard-like interaction term Hu = 
Und^fidi only to a single d-orbital of the Go 3c?-shell 
where Uda- is the number operator for the d-orbital and a 
spin a. Since the Goulomb interaction in the correlated 
d-orbital has already been taken into account on a mean- 
field level in the KS-DFT calculation, a double-counting 
correction (DGG) term has to be subtracted from the KS 
Hamiltonian projected onto the d-orbital = {d\ |d): 

= ^d ~ ^dc- (2) 

In contrast to previous work the DGG is chosen such 
that a certain occupancy is achieved, i.e. for achieving 
particle-hole (ph) symmetry {ud = 1) we choose Cdc such 
that Cd = —Ujl. Note that ph symmetry is only ap¬ 
proximately achieved since the coupling of the d-orbital 
to the rest of the system (see below) is generally not ph 
symmetric. 

The interacting d-orbital coupled to the electronic bath 
given by the rest of the system (i.e. substrate and tip) de¬ 
fines an Anderson impurity model (AIM)I^. An effective 
description of the coupling of the d-orbital to the bath 
is given by the so-called hybridization function A.d{u}) 
which can be obtained from the KS GF by 


We consider a magnetic atom (here: Go) that is placed 
on a metallic substrate (here: Gu(OOI) or Gu(IlI)). A 
Gu STM tip is placed directly above the Go atom gA 
away so that we are in the tunneling regime. The system 
is divided into three parts as shown in Fig. two metal 
leads S and T, representing the bulk electrodes connected 
to the substrate and STM tip, respectively, and the de¬ 
vice region (D) which contains the magnetic atom and 
part of the surface and the STM tip. 

We perform DFT based ab initio quantum transport 
calculations using the ANT.G packag^^ The electronic 
structure of the D region is calculated on the level of 
Kohn-Sham (KS) DFT emplo ying t he LSDA functionaP^ 
in the SVWN parametrizatiorPZEH and a minimal Gaus¬ 
sian basis set including the valence (4s4p3d) and outer 


Ad(uj) — uj + g, - 6° - [Gd{uj)] , (3) 

where is the KS GF projected onto the d-orbital, i.e. 
G^ = (d| Gd M)- The imaginary part of Ad yields the 
broadening F^; of the d-orbital due to the coupling to the 
rest of the system. 

The AIM is now solved in the OGA.I^ The solution 
yields the self-energy E(i(w) describing the dynamic cor¬ 
relations of the d-orbital. The correlated GF of the d- 
orbital is then given by 

Gd{u}) = ([G2(w)] ^ — T,d{uj) + Cdc) • (4) 

Its imaginary part yields the spectral function or LDOS 
of the d-orbital Pdi^) = —IniGd(a;)/7r. Gorrespondingly, 
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we obtain the correlated GF for the D region as: 


Gd(cc) = - (S^ - edc)Pd 


( 5 ) 


This allows us to calculate the transmission function us¬ 
ing the Caroli expressiorP^, 


T{lo) = Tr 


GofTCtfs 


( 6 ) 


where the coupling matrices for the leads are defined by 
fx/S = * (^T/S - ^T/s) ■ 

The self energies Sx/s are typically symmetric, so that 
the coupling matrices are twice the imaginary part of 
the self energies. For low temperature and small bias 
voltages, current and conductance can be related t o the 
transmission function using the Landauer formulaP^^. 
For the typical STM setup considered here most of the 
applied bias voltage drops at the STM tip. In that case 
the conductance is simply given by 

G{V) = ‘^neV). (8) 

We note that the use of the Landauer formula for the 
conductance is justified in the limit of small bias volt¬ 
ages compared to the Kondo temperature. In this limit 
transport occurs via the Kondo resonance and thus is es¬ 
sentially one-body like (apart from renormalization) and 
phase coherent so that the full non-equilibrium expres¬ 
sion for transport through an interacting region given W 
Meir-Wingreen reduces to the simpler Landauer result.!^ 
For larger bias voltages, deviations from the Landauer re¬ 
sult can occuil^^^^, and one would have to make use of 
the Meir-Wingreen equatioiP^ which requires the solu¬ 
tion of the AIM out of equilibrium. 


III. FANO-LORENTZ AND FANO-FROTA LINE 
SHAPES 

Fano line shapes or resonances, originally introduced 
by Fano in the context of autoionization and elastic elec¬ 
tron scattering by heliunP^, generally arise in resonant 
scattering processes due to quantum interference between 
a quasi discrete resonant state and a broad background 
continuum. The interference leads to an asymmetric line 
shape in the scattering cross section at energies close to 
the resonance energy that is well described by the Fano 
function 


/(e) oc 


(g + e)^ 

£2 + 1 ’ 


( 9 ) 


where the parameter q controls the shape of the Fano 
function, and e is the energy with respect to the reso¬ 
nant level. Eq. [^can also be obtained from the complex 


representation of a Lorentzian multiplied by a phase fac¬ 
tor 


Pfl(w) = Im 



A 


Ct) — Wo + IF 


+ Po, 


( 10 ) 


where A is the amplitude, F is the half-width of the 
Lorentzian, uq the resonance energy and po a constant 
offset. Using q = tan(0g/2), e = (w — wo)/F and some 
algebra, this Fano-Lorentz (FL) line shape can be shown 
to be equivalent to the original Fano formula (see the 
Appendix): 


Pfl(w) 


A( {q + er \ 
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STM spectroscopy of Kondo impurities presents a sim¬ 
ilar situation: The STM tip probes the continuous con¬ 
duction electron density of states which interacts with 
the Kondo resonance at the Fermi level. The interfer¬ 
ence of different tunneling paths then leads to Fano- 
type line shapes in the conductance spectra. Assuming a 
Lorentzian form for the Kondo resonance naturally leads 
to Fano-Lorentz line shapes given by ( [Tl] ). However, in 
Refs, [m and Frota showed that the Kondo peak is 
actually better described by a line shape now known as 
a Frota line shape: 


/^Frota(^) — A • Re 


iF p 

U! — UlQ + F 


( 12 ) 


where the Frota parameter Ff is related to the actual 
half-width F of the resonance by F ^ 2.54 Ff. A is the 
amplitude and ujq the position of the Frota resonance. In 


analogy with Eq. (10) we dehne a Eano-Erota (FF) line 
shape as a generalized Frota curv^^H for describing the 
transmission function close to the Kondo resonance 


Tff(w) = —A ■ Re 




zFf 


w — Wo -f IFf 


+ To. (13) 


The phase (j)q has the same meaning in the Lorentz and 
in the Frota case: A value oi (j)q = Q leads to a dip, 
= TT to a peak and (()q = 7r/2 to a symmetric Fano 
line shape. In Fig. we compare Fano-Frota and Fano- 
Lorentz features, choosing identical amplitudes and half¬ 
widths. Note that for the same half-width Frota line 
shapes have a slower decay than the Lorentzian ones. 

In Fig. we show FL and FF fits to the Kondo peak 
(left) in the calculated spectral function pd(oj) and the 
corresponding Fano line shape (right) in the calculated 
transmission function T(w) for the case of the z^-orbital 
for the Co on Cu(OOl) system, discussed in detail in the 
following section. For both spectral and transmission 
function, the resonance center is well-described by FF 
and FL fits. However, only the FF ht yields an accurate 
description of the flanks and the long range decay. In the 
following, we will therefore use Eqn.[^to fit transmission 
functions. 
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FIG. 2. Comparison of Fano-Frota and Fano-Lorentz line 
shapes for different values of the q parameter but for iden¬ 
tical amplitudes and half-widths. 
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FIG. 3. Fano-Lorentz (FL) and Fano-Frota (FF) fits of the 
impurity spectral function pd(tu) (left) and the transmission 
function T{uj) (right) for a^-orbital, Co@Cu(001), U = 2 eV 
and td = —1.0 eV. See Sec. IIV A] for details. 




IV. RESULTS 

A. Co adatom on Cu(OOl) surface 

The system under consideration is shown in the left 
panel of Fig|^ A Cobalt atom is deposited at the hollow 
site of a Cu(OOl) surface. The Cu(OOl) surface is modeled 
by three Cu slabs of 36, 25 and 16 atoms, respectively, 
which are embedded into a Bethe lattice to describe the 
inhnitely extended surface. We model the STM tip by a 
small pyramid of Cu atoms grown in the (001) direction, 
also embedded into a Bethe lattice. The tip is placed 
directly above the Co atom in a distance of 6A, so that 
the system is in the tunneling regime. 

As explained in Sec. [^we now compute the hybridiza¬ 
tion functions of the Co 3d-orbitals (see right panel of 
Fig. 1^. The four-fold symmetry of the Cu(OOl) sur¬ 
face leads to a splitting into four groups. The xz- and 
?/ 2 ;-orbitals are degenerate (in the following, results for 
the j/z-orbital are omitted) and exhibit the strongest hy¬ 
bridization at the Fermi level. The hybridization func¬ 
tions of and have comparable values around the 
Fermi level. The xy-orbital has the lowest hybridization 
in the displayed energy window. All hybridization func¬ 
tions show a moderate energy dependence. Note that the 
hopping between different Co 3d orbitals is zero, i.e. they 
do not couple to each other on the single-particle level. 

While the hybridization function is calculated ab ini¬ 



FIG. 4. Left: Geometry of the Co atom deposited on a 
Cu(OOl) surface; dark gray/red: Co, gray/yellow: Cu, light 
gray/blue: Bethe lattice. Right: Imaginary part of the hy¬ 
bridization function for the Co 3d shell. 


tio, the Coulomb interaction U is used as a parameter 
that allows us to tune the Kondo coupling strength and 
explore the effect of the width of the Kondo peak on the 
transmission line shape. But in order to have an estimate 
of the magnitude, we have also calculated U ab initio for 
each of the d-orbitals by constrained RPA calculations as 
described in Ref.^^. We find values for U ranging from 
1.8 eV to 2.6 eV.I^ Accordingly, we choose the U param¬ 
eters to vary between 2 eV and 3 eV. 

The hybridization functions from Fig. [^together with 
the energy level and the effective Coulomb interac¬ 
tion U define an AIM which is solved in the OCA^^. It 
is a known issue of OCA that at too low temperatures 
(1-2 orders of magnitude below Tk) it gives rise to spu¬ 
rious non-Fermi liquid behavior and related artifacts in 
the impurity spectral function, leading to an overestima¬ 
tion of the height of the Kondo peak and an unphysical 
self-energy with positive imaginary parlP^. We circum¬ 
vent this problem by lowering the temperature only to 
the point where the imaginary part of the self-energy be¬ 
comes zero. At this point Fermi liquid behavior is obeyed, 
and the unitary limit of the Kondo peak is exactly recov¬ 
ered. 

Fig. shows impurity spectral functions Pd{oj) of the 
z^-orbital for different values of the AIM parameters 
and U. For = —U /2 (red solid and blue dotted curves) 
we have approximate particle-hole symmetry: the Kondo 
peak is centered close to, but slightly above the Fermi 
level. Note that exact particle-hole symmetry is not 
achieved because of the non-constant hybridization func¬ 
tion. As expected, when U is increased the Kondo tem¬ 
perature and hence the width of the Kondo peak decrease 
strongly. On the other hand detuning the system from 
particle-hole symmetry by shifting ea leads to a strong 
increase of the Kondo temperature due to charge fluctu¬ 
ations (green dashed, magenta dashed-dotted curves). 

We now calculate the correlated transmission func¬ 
tions for Kondo peaks in different d-orbitals. Fig. 
shows transmission line shapes for different d-orbitals for 
U = 2 eV and = —1.0 eV. In order to make the fea¬ 
tures more clearly visible, here and in the following the 
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FIG. 5. Impurity spectral functions for the orbital of 
Co@Cu(001) for different Anderson impurity model param¬ 
eters U, ed. 


transmission background was subtracted.^ We find that 
the line shapes are indeed different for each orbital. We 
observe approximately antisymmetric Fano line shapes 
(g « 1) for and — y^, and more peak-like feature 
(g ^ 1) for xz and xy. In order to quantitatively describe 
the line shapes, we perform Frota fits to determine the 
g parameter and width of the line shapes, as explained 
before in Sec. Ill The and x^ — y^ orbitals have com¬ 
parable g values of 1.52 and 1.72, respectively. For xz, 
q becomes negative (—3.9) and for xy we find the most 
pronounced peak with g = 4.95. The width of the Fano 
features differs significantly, and in accordance with their 
hybridization strength at the Fermi level. Note that a 
feature with a very small width, as e.g., in the case of 
xy, might never be observed in an actual experiment, 
because of the Kondo temperature being much too low 
and because of limited resolution. 


We now vary the Coulomb repulsion U and introduce 
charge fluctuations by shifting the d-level position e^, as 
can be seen in Fig. When varying U, but maintaining 
particle-hole symmetry, the actual shape of the transmis¬ 
sion features is only weakly affected, while the widths of 
the features change strongly, as has already been seen 
and discussed for the spectral functions in Fig.[^ When 
introducing charge fluctuations, the Kondo peak becomes 
asymmetric (see Fig.[^. This asymmetry is also reflected 
in the transmission line shapes. We find that the g pa¬ 
rameter consistently increases when is shifted down¬ 
wards. For positive g (z^, — 2/^, xy) lowering makes 

the line shapes more peak-like, while for negative g (xz), 
lowering ed leads to more dip-like line shapes. 

Hence while the choice of AIM parameters U and Cd 
does affect the transmission line shapes to some degree, it 
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FIG. 6. Transmission functions for different d-orbitals of 
Go@Gu(001). Goulomb repulsion U = 2 eV, and energy level 
= — 1 eV (approximate particle-hole symmetry). The red 
continuous curves show the calculated transmission, the black 
dashed curves Fano-Frota fits. The transmission background 
has been subtracted.^ 


does not completely change its symmetry. For example, 
the sign of the g factor does not change. 

While the signal width is determined by the hybridiza¬ 
tion and choice of AIM parameters exclusively, the signal 
amplitude decisively depends on the system geometry. 
Because we chose the z-axis as our transport direction, 
a Kondo peak in the z^-orbital results in a much more 
dominant feature compared to the remaining d-orbitals, 
as can be seen in Figs. and Hence if there is a 
Kondo peak in the z^-orbital, the corresponding Fano 
feature will dominate in the transmission regardless of 
what happens in the other orbitals. Also Fano features 
due to Kondo peaks in orbitals other than the z^-orbital 
might be difficult to discern from the background if the 
background dispersion is strong compared to the Fano 
amplitudes. This statement remains true even if the STM 
tip is shifted laterally by moderate distances of a few A. 
Although tunneling into orbitals other than z^ becomes 
more favorable upon a lateral shift of the tip, the fea¬ 
ture due to the Kondo peak in the z^ remains the most 
dominant one. 


B. Co adatom on Cu(lll) surface 

The next system we focus on is a Cobalt atom, de¬ 
posited at the ’hep’ hollow site of a Cu(lll) surface, as 
can be seen in the left panel of Fig. The surface is 
modeled by three Cu slabs of 27, 37 and 27 atoms, re¬ 
spectively, which are connected to a Bethe lattice. The 
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FIG. 7. Transmission functions for different d-orbitals of 
Co@Cn(001). Coulomb repulsion U = 3 eV, vary occupation 
by shifting ed. The red continuous curves show the calculated 
transmission, the black dashed curves Fano-Frota fits. The 
transmission background has been subtracted.^ 





FIG. 8. Left: Geometry of the Co atom deposited on a 
Cu(lll) surface; dark gray/red: Co, gray/yellow: Cu, light 
gray/blue: Bethe lattice. . Right: Imaginary part of the 
hybridization function of the Co Sd-shell. 




FIG. 9. Transmission functions for different d-orbitals of 
Co@Cu(lll). Coulomb repulsion U = 2 eV, td = —1.0 eV. 
The red continuous curves show the calculated transmission, 
the black dashed curves Fano-Frota fits. The transmission 
background has been subtracted.^ 


tip is described by a Cu(lll) pyramid, consisting of 10 
copper atoms, also connected to a Bethe lattice. The 
threefold symmetry splits the hve orbitals of the Co 3d- 
shell into three groups: the non-degenerate z^-orbital 
(to = 0) and two doubly degenerate groups, one with 
TO = ±1 {xz- and j/z-orbitals) and one with to = ±2 
{xy- and x'^ — j/^-orbitals). The right panel of Fig. 
shows the hybridization functions for each of the three 
groups. The group with the xz- and yz-orbitals exhibit 
the strongest hybridization at the Fermi level, the group 
with the x'^ — y'^- and xy-orbitals the weakest. 

We proceed as described in the previous section and 
calculate transmission functions for the d-orbitals of 
Co@Cu(lll), assuming a Coulomb repulsion oiU = 2 eV 
and (approximate) particle-hole-symmetry ed = —1.0 eV 
(Fig. §. Again, we find different line shapes for each 
orbital. The z^ orbital gives the most peak-like trans¬ 
mission feature with q = 2.77, for — y^ we observe a 
transmission peak with a negative q = —2.44. The xz- 
orbital results in a Fano-type feature with q = —1.61. 
The widths of the transmission features differ consid¬ 
erably, with the xz- and yz-orbitals having the largest 


width, and the xy- and — y^-orbitals the lowest. The 
z^-orbital again has the highest signal amplitude, as it 
couples strongly to the tip conduction electrons. 


In Fig. we calculate line shapes for different AIM 
parameters U and Cd- We observe a similar behavior as 
for Co@Cu(001). When staying in the particle-hole sym¬ 
metric case and increasing U (middle column of Fig. 101, 
the line shapes remain similar, with slightly increased q 
values. We introduce charge fluctuations by shifting the 
position of ed (left and right column of Fig. 10). The q 
parameter increases when moving to lower energies. 
For positive q values, as for this leads to more peak¬ 
like line shapes, while for negative q values, as for xz and 
x^ — y^, it leads to more fano- or dip-like line shapes. The 
only exception to this behavior occurs for the xz orbital, 
U = 3 eV and ed = —0.8 eV. It has a very high Kondo 
temperature and equivalently wide Fano feature, and the 
Fano-Frota fit fails for negative energies. This suggests 
that the Fano line shape overlaps with other transmission 
features that alter the final line shape. 
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FIG. 10. Transmission functions for different d-orbitals of 
Co@Cu(lll). Coulomb repulsion [7 = 3 eV, vary occupation 
by shifting td- The red continuous curves show the calculated 
transmission, the black dashed curves Fano-Frota fits. The 
transmission background has been subtracted.!^ 
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FIG. 11. Temperature dependence of two different line shapes 
for Co@Cu(lll), xz, U = 3 eV, ed = —1.5 eV and = 
— 2.2 eV, respectively. Top: Transmission. Bottom: q param¬ 
eter; the lines are a guide for the eye. 


C. Temperature dependence and the occurrence of 

dips 


The results presented so far are for the case of T 


0 (according to the criterion discussed in Sec. IVA). 


We now study the temperature dependence of two line 
shapes: One tending towards a peak {q > 1) and one 
tending towards a dip {q < 1). We pick the xz orbital 
of Co@Cu(lll), U = 3 eV, €d = -1.5 eV {q = -1.49) 
and Cd = —2.2 eV {q = —0.87), respectively. The top 
row of Fig. shows the evolution of the aforementioned 
two line shapes. For increasing temperature, the signal 
amplitude diminishes, while its width grows. The peak 
does not decay symmetrically. The ’peak’ component 
of the Fano feature decays faster than the ’dip’ compo¬ 
nent of the feature, so that, in both cases, the feature 
as a whole becomes increasingly dip-like with increas¬ 
ing temperature. In order to quantify that, we perform 
Fano-Frota fits and calculate the q parameter. We find 
that the q parameter decreases considerably when tem¬ 
perature is rising, irrespective if the feature tends more 
towards peak or dip in the T —>• 0 case. 


V. A SIMPLIFIED MODEL 

The interference mechanism leading to d ifferent Fano 
line shapes still is a matter of discussioiP^J^Iinill], 
expand on this discussion by introducing a simple model 


that allows us to determine transmission line shapes from 
ab initio parameters. Fig. |12| shows a schematic drawing 
of our model system. The central assumption is that the 
quantum interference primarily occurs on the magnetic 
adatom, namely between one s-type and/or p-type level 
(in the following, we will simply call it the conduction 
level c) and the correlated d-level. Both levels are in 
contact to the tip T and the surface S, and the respective 
interactions are taken into account by coupling matrices 
Fx/s- As a second central assumption we neglect the 
direct tunneling from the tip to the surface. 

The starting point of our model is the correlated 
Green’s function of the effective atom comprising the con¬ 
duction c-level and the correlated d-level of the magnetic 
atom. 

Gk(.uj) = (wPa - Pa - A(w) - Sd(a;)Pd)”' 

w - Cc - Ac(w) -Vcd - 
-Vdc - ^dc{uj) ui -ed - Ad(a;) - ^^(a;) 

(14) 

Pa is a projector onto the effective atom A, while Pd 
projects onto the d-level only. All parameters can either 
be extracted from the KS-calculation (e^, Cc, 14d, A(w)) 
or from the OCA-calculation (I](i(a;)), while the chemi¬ 
cal potential has been set to zero /i = 0. The diagonal 
elements of the hybridization function A(a;) lead to a 
shift (real part) of the level position of Ec and e^, respec- 




























tip 



FIG. 12. Sketch of the simplified model. The effective atom 
A is described by the correlated d-level and one conduction 
electron level c, in contact with the surface and the tip. 


lively, and yield an effective level broadening (imaginary 
part). Also note that the hybridization function has off- 
diagonal components Acd{ui) = Adc(oj), which can be 
understood as an additional hopping between the c- and 
d-level mediated by hoppings via the substrate, to give a 
total effective coupling of Vcd = Vcd + Acd- The coupling 
matrices rg/x(w) necessary for calculating the transmis¬ 
sion function by Eq. can be obtained by decomposing 
the hybridization function into a tip (At(w)) and a sur¬ 
face (As(w)) component and taking the imaginary parts, 
be. rs/T(‘^) = —2ImAs/T(w). 

For the conduction level c of the effective atom we 
choose the s- or p-orbital that couples to the correlated 
d-orbital. In the case of the (001) and the (111) sub¬ 
strates the z^-orbital couples to the s- as well as the 
P 2 -orbital. In this case we apply a unitary transforma¬ 
tion in the subspace of the s- and Pz-orbitals such that 
the z^-orbital decouples completely from one of the or¬ 
bitals in the new basis. The sp^-hybridized orbital cou¬ 
pling to the is then found to be the linear combination 
\spz) oc Vsz 2 |s) -b \pz) where Vgz'^ and Vp^z'^ are the 
effective hoppings of the z^-orbital with the s- and Pz- 
orbitals, respectively. On both surfaces, the a;z-orbital 
couples to Px and the yz-orbital to Py. For the (001) 
surface both the and the xy-orbitals do not in¬ 

teract with any of the s- or p-orbitals on the atom, while 
on the (111) surface, they do interact with the Py- and 
Pa;-orbitals, respectively. 


In Fig. [T^and[l4l we compare line shapes calculated 
for the simplified model with the full ab initio results 
from Sec. IV For the Co@Cu(001) surface, the simplified 
model consisting of the z^-orbital and the sp 2 -hybridized 
orbital reproduces the line shape of the z^-orbital quite 
well. Only the peak character is slightly overestimated. 
In the case of the ccz-orbital the line shape of the simpli¬ 
fied model including the Pa;-orbital is in excellent agree¬ 
ment with that of the full ab initio calculation. For 
the — y^-orbital the agreement between the simpli¬ 
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FIG. 13. Transmissions calculated ah initio with the 
ANT.G package (see Sec. |IV[ ) and for the simplified model; 
Go@Gu(001), [7 = 2 eV, ed = — 1 eV. The transmission func¬ 
tions are rescaled and offset for better visibility. 


fied model and the full calculation is not as good. As 
stated before this orbital does not interact with any s- 
or p-orbital on the Co atom. Hence the transmission of 
the simplified model reproduces simply the Kondo peak 
in the spectral function since no interference is taking 
place. On the other hand the full transmission shows a 
somewhat asymmetric Fano feature {q « 1.7) indicating 
that interference with some substrate state(s) must take 
place, which is not included in the model. Finally, for the 
ccy-orbital we find very good agreement between the sim¬ 
plified model and the full calculation. The line shape in 
both calculations simply reproduces the Kondo peak in 
the spectral function of the a:?/-orbital indicating the ab¬ 
sence of any interference effects between this d-level and 
s- and p-levels on the atoms as well as substrate states. 

We find a somewhat similar picture for Co@Cu(lll). 
For the xz-orbital the model including the interaction 
with the p 2 ,-orbital gives a line shape in excellent agree¬ 
ment with the full calculation. Also for the the 

simplified model including the p^-orbital on the atom re¬ 
produces the line shape of the full calculation very well. 
However, in the case of the z^-orbital the simplified model 
including the spj-hybridized orbital fails quite badly in 
reproducing the line shape of the full calculation. Ap¬ 
parently, interference with tunneling paths to substrate 
states play an important role here. 
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FIG. 14. Transmissions calculated ab initio with the 
ANT.G package (see Sec. IVI and for the simplified model; 
Go@Gu(lll), U = 2 eV, Cd = — 1 eV. The transmission func¬ 
tions are rescaled and offset for better visibility. 


VI. DISCUSSION 


For Co@Cu(001), we found transmission line shapes 
ranging from asymmetric Fano features with positive , 
— y^) and negative (xz) q values to a more peak¬ 
like feature (xy). The line shapes are determined by 
the interference of different tunneling paths. Our sim¬ 
plified model calculations indicate that for and xz 
the interference takes place on the adatom between the 
correlated d-level and the non-interacting sp-levels cou¬ 
pling to the d-orbital. For the xj/-orbital, no interference 
occurs between the conduction and impurity tunneling 
channels. Hence one directly observes the shape of the 
Kondo peak in the transmission. On the other hand, for 
the x'^ — y^-orbital, the interference mechanism proba¬ 
bly involves the Cu substrate states which are not cap¬ 
tured by the simplified model. Experimentally, asym¬ 
metric Fano line shapes were reported with q ^ 1.1 — 1.2 
in the tunneling regirnd^^^^^. The measured line shapes 
are comparable to the features we found both in the 
z^- and x'^ — y^-orbitals (see Figs. and [^, although 
the z^-orbital yields a slightly better agreement. Bet¬ 
ter agreement with experiment can surely be achieved 
by adjusting the Anderson model parameters and fitting 
the calculated spectra with the experimental ones. We 
would like to stress though that finding good agreement 
with experiment is not the primary goal of this work, but 
rather to demonstrate how different orbital symmetries 
give rise to different Fano-Kondo line shapes. A recent 
study by one of u^^ found an underscreened Kondo ef¬ 
fect for Co@Cu(001), where the z^ and x^ — y^ are nearly 
half filled, but only the z^ orbital is Kondo screened at fi¬ 
nite temperatures due to its higher Kondo temperature. 
Ref. 131] comes to similar conclusions, finding a Kondo 
peak in the z^-orbital with q = 1.2 in the tunneling 
regime and explaining it due to the interference of the 
z^- with the s-orbital. 

For Co@Cu(lll), we found asymmetric to peak-like 
Fano line shapes with positive (z^) and negative (xz, x^ — 
y'^) q values. For the latter two, we can understand the 
tunneling interference in terms of the model presented 
in the previous section. The interference occurs on the 


magnetic atom, between the conduction electron channel, 
modeled by one of the p-orbitals, and the respective d- 
level. For z^, which is interacting with the sp^ hybridized 
level, our model fails, indicating that interference with 
substrate states plays an important role here. 

Exp eriment ally, dips were reported with q values close 
to zerd ^^bz lM] does not seem to agree with any of 

the calculated line shapes. The z^-orbital, aligned in the 
transport direction, again shows the strongest signal, but 
is rather peak-like. The closest candidate to a dip-like 
line shape is the a:z-orbital, particularly when increas¬ 
ing the occupancy relative to half-filling by moving the 
d-level position downwards in energy (see Fig. 10). In 
Sec. IV C we studied its temperature dependence, and 
found that the line shape became increasingly dip-like 
when increasing temperature. However, note that in our 
calculations for the xz-orbitals we find q < 0 while in 
experiment q is always positive. 

Probably, the surface state of the Cu(lll) surfacd^ 
plays an imp ortant role for determining the line 
shapdSMMlllI since its tunneling amplitude may be 
twice as strong compared to tunneling into bulk stated^. 
However, our embedded cluster calculation probably does 
not capture the surface state properly. The importance of 
the surface state for reproducing the correct line shape in 
the Co@Cu(lll) system is also stressed in Ref. ED where 
the surface state is not properly captured and the correct 
q value could not be reproduced either. 


VII. CONCLUSIONS 

In summary, we have calculated the orbital signatures 
of Kondo peaks in the STM spectra of transition metal 
adatom systems, namely Co@Cu(001) and Co@Cu(lll). 
Our calculations show that the measured line shapes al¬ 
low us to draw some conclusions on the d-orbital(s) in¬ 
volved in the Kondo effect since the line shape depends 
to a large extent on the coupling of the d-orbital to the 
sp-orbitals on the adatom, which in turn is determined 
by the orbital symmetry. However, also temperature, ef¬ 
fective interaction (7 and in particular the occupancy of 
the d-orbital have a strong influence on the actual line 
shapes. Also, if multi-orbital effects are important for 
the actual shape of a Fano-Kondo feature, this approach 
per se is not appropriate. Nevertheless, even in the case 
of a multi-orbital Kondo effect, often one orbital will be 
dominant in the tunnel spectra. In fact if a Kondo res¬ 
onance forms in the z^-orbital, the corresponding Fano 
feature will be dominant in the tunnel spectrum for the 
typical case of an s-type STM tip, so that Kondo features 
coming from other d-orbitals are likely not visible. These 
results are also relevant for STS of transition metal com¬ 
plexes on metallic substrate^^^^^, maybe even more so 
since tunneling into surface states is less important there. 

We stress that the here developed method can in prin¬ 
ciple also be applied to the contact regime. However, 
unlike in the tunneling case, in the contact regime the 
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voltage can no longer be assumed to mainly drop be¬ 
tween tip and adatom. Rather, the voltage drop will 
distribute in some way over the contact according to the 
actual geometry of the contact regioiJ^, and needs to be 
calculated or estimated. Moreover, the actual contact ge¬ 
ometry is probably also relevant for the coupling between 
d-orbitals and conduction electrons and thus also has a 
strong influence on the line shapes. Therefore possible 
contact geometries need to be explored and relaxed with 
some care. 

Based on our results, we propose a poor man’s method 
to obtain information on the orbital(s) involved in the 
Kondo effect measured in an actual experiment solely on 
the basis of a density functional theory calculation of the 
system: by tailoring an appropriate self-energy for each 
orbital such that the width of the resulting Kondo peak in 
that orbital reproduces the width of the measured Fano- 
Kondo line shape, one can calculate the corresponding 
line shapes and compare to experiment. 
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Appendix A: Complex and real Fano line shapes 


Here we derive the real Fano line shape given by 
Eq. (11) from its complex representation in Eq. (10): 


Pfl(w) = Im 




A 


^uj — ujQ + ir 

Introducing the abbreviation e = (co — a;o)/F, we have 

e — i 


(Al) 


Pfl = p • Im 


(cos((/)g) +isin(0g)) ^ 


A 

r 

A 

f 

A 

f 


-cos{(j)q) +esm{(j)q) 


— cos^(^) -I- sin^(^) -I- 2esin(^) cos(^) 
-|- 1 


-1 -I- tan^(^) -I- 2etan(^) 


-I- 1 


cos 


Defining q = tan((}()g/2), we arrive at 


A 


Pfl = 


(9 + e)^ 


- 1 


1 + q^ 


(A2) 


which is the same as Eq. 


A. Hewson, The Kondo Problem to Heavy Fermions, edited 
by D. Edwards and D. Melville (Cambridge University 
Press, Cambridge, 1997). 

^ W. de Haas, J. de Boer, and G. van dn Berg, Physica 1 , 
1115 (1934) 

^ M. P. Sarachik, E. Corenzwit, and L. D. Longinotti, Phys. 
Rev. 135, A1041 (1964) 

J. Kondo, Progress of Theoretical Physics 32, 37 (1964) 

® A. Posazhennikova, B. Bayani, and P. Coleman, Phys. 
Rev. B 75, 245329 (2007) 

° N. Roch, S. Florens, V. Bouchiat, W. Wernsdorfer, and 
F. Balestro, Nature 453, 633 (2008) 

^ J. Nygard, D. H. Cobden, and P. E. Lindelof, Nature 408, 
342 (2000) 

“ P. Jarillo-Herrero, J. Kong, H. S. van der Zant, C. Dekker, 

L. P. Kouwenhoven, and S. De Franceschi, Nature 434, 
484 (2005) 

“ J. Li, W.-D. Schneider, R. Berndt, and B. Delley, Phys. 
Rev. Lett. 80, 2893 (1998) 

V. Madhavan, W. Chen, T. Jamneala, M. F. Crommie, 
and N. S. Wingreen, Science 280, 567 (1998) 

H. C. Manoharan, C. P. Lutz, and D. M. Eigler, Nature 
403, 512 (2000) 

N. Knorr, M. A. Schneider, L. Diekhoner, P. Wahl, and 

K. Kern, Phys. Rev. Lett. 88, 096804 (2002) 

K. Nagaoka, T. Jamneala, M. Grobis, and M. F. Crommie, 
Phys. Rev. Lett. 88, 077205 (2002) 

A. J. Heinrich, J. A. Gupta, G. P. Lutz, and D. M. Eigler, 
Science 306, 466 (2004) 


A. Zhao, Q. Li, L. Ghen, H. Xiang, W. Wang, S. Pan, 

B. Wang, X. Xiao, J. Yang, J. G. Hou, and Q. Zhu, Science 
309, 1542 (2005) 

M. F. Crommie, Science 309, 1501 (2005) 

V. lancu, A. Deshpande, and S.-W. Hla, Phys. Rev. Lett. 
97, 266603 (2006) 

N. Neel, J. Kroger, L. Limot, K. Palotas, W. A. Hofer, 
and R. Berndt, Phys. Rev. Lett. 98, 016801 (2007) 

U. Fano, Phys. Rev. 124, 1866 (1961) 

H. Priiser, M. Wenderoth, A. Weismann, and R. G. Ul- 
brich, Phys. Rev. Lett. 108, 166604 (2012) 

H. O. Frota and L. N. Oliveira, Phys. Rev. B 33, 7871 
(1986) 

H. O. Frota, Phys. Rev. B 45, 1096 (1992) 

M. Plihal and J. W. Gadzuk, Phys. Rev. B 63, 085404 

( 2001 ) 

V. Madhavan, W. Ghen, T. Jamneala, M. F. Grommie, 
and N. S. Wingreen, Phys. Rev. B 64, 165412 (2001) 

C. -Y. Lin, A. H. Castro Neto, and B. A. Jones, Phys. Rev. 
Lett. 97, 156102 (2006) 

O. IJjsaghy, J. Kroha, L. Szunyogh, and A. Zawadowski, 
Phys. Rev. Lett. 85, 2557 (2000) 

A. Schiller and S. Hershfield, Phys. Rev. B 61, 9036 (2000) 

P. Wahl, L. Diekhoner, M. A. Schneider, L. Vitali, G. Wit- 
tich, and K. Kern, Phys. Rev. Lett. 93, 176603 (2004) 

J. Merino and O. Gunnarsson, Phys. Rev. B 69, 115404 
(2004) 

J. Merino and O. Gunnarsson, Phys. Rev. Lett. 93, 156601 
(2004) 






























11 


P. P. Baruselli, R. Requist, A. Smogunov, M. Fabrizio, and 
E. Tosatti, Phys. Rev. B 92, 045119 (2015) 

L. Vitali, R. Ohmann, S. Stepanow, P. Gambardella, 
K. Tao, R. Huang, V. S. Stepanyuk, P. Bruno, and 
K. Kern, Phys. Rev. Lett. 101, 216802 (2008) 

N. Neel, J. Kroger, and R. Berndt, Phys. Rev. B 82, 
233401 (2010) 

D.-J. Choi, M. V. Rastei, P. Simon, and L. Limot, Phys. 


Rev. Lett. 108, 266803 (2012) 


D. Jacob and J. J. Palacios, 

The Journal of Chemical 


Physics 134, 044118 (2011) 

W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965) 

C. J. Slater, The Self-Consistent Field for Molecular and 
Solids, Quantum Theory of Molecular and Solids, Vol. 4 
(McGraw-Hill, New York, 1974). 

S. H. Vosko, L. Wilk, and M. Nusair, Canadian Journal 
of Physics 58, 1200 (1980), 

W. J. Hehre, R. F. Stewart, and J. A. Pople, The Journal 
of Chemical Physics 51, 2657 (1969) 

J. B. Collins, P. von R. Schleyer, J. S. Binkley, and J. A. 
Pople, The Journal of Chemical Physics 64, 5142 (1976) 
P. J. Hay and W. R. Wadt, The Journal of Chemical 
Physics 82, 270 (1985) 

W. R. Wadt and P. J. Hay, The Journal of Chemical 
Physics 82, 284 (1985) 

P. J. Hay and W. R. Wadt, The Journal of Chemical 
Physics 82, 299 (1985) 

H. A. Bethe, Proc. Roy. Soc. London A: Math. Phys. Eng. 
Sci. 150, 552 (1935) 

K. Haule, S. Kirchner, J. Kroha, and P. Wolfle, Phys. Rev. 
B 64, 155111 (2001) 

^ D. Jacob, K. Haule, and G. Kotliar, Phys. Rev. Lett. 103, 
016803 (2009) 

D. Jacob, Journal of Physics: Condensed Matter 27, 
245606 (2015) 

P. W. Anderson, Phys. Rev. 124, 41 (1961) 

C. Caroli, R. Combescot, D. Lederer, P. Nozieres, and 

D. Saint-James, Journal of Physics C: Solid State Physics 
4, 2598 (1971) 


R. Landauer, IBM Journal of Research and Development 
1, 223 (1957) 

S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge University Press, Cambridge, 1995) Cam¬ 
bridge Books Online. 

Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992) 


M. H. Hettler, J. Kroha, and S. Hershfield, Phys. Rev. B 
58, 5649 (1998). 


C. A. Balseiro, G. Usaj, and M. J. Sanchez, 

Journal of 

Physics: Condensed Matter 22, 425602 (2010) 



P. Roura-Bas, Phys. Rev. B 81, 155327 (2010). 


H. Priiser, M. Wenderoth, P. E. Dargel, A. Weismann, 
R. Peters, T. Pruschke, and R. G. Ulbrich, Nat Phys 7, 
203 (2011) 

More specifically, we obtain 1.80 eV for the 2 ;^-orbital, 

I. 78 eV for the xz- and j/ 2 -orbitals, 1.81 eV for the 
orbital and 2.59 eV for the xy orbital. 

T. A. Costi, J. Kroha, and P. Wolfle, Phys. Rev. B 53, 
1850 (1996) 

The background is calculated by calculating the transmis¬ 
sion function without adding the self-energy but pushing 
the respective d level away from the Fermi level. 

H. G. Luo, T. Xiang, X. Q. Wang, Z. B. Su, and L. Yu, 
Phys. Rev. Lett. 92, 256602 (2004) 

R. Zitko and T. Pruschke, New Journal of Physics 12, 
063040 (2010) 

M. F. Crommie, C. P. Lutz, and D. M. Eigler, Nature 
363, 524 (1993) 

A. Bogicevic, S. Ovesson, P. Hyldgaard, B. 1. Lundqvist, 
H. Brune, and D. R. Jennison, Phys. Rev. Lett. 85, 1910 
( 2000 ) 

O. Jeandupeux, L. Biirgi, A. Hirstein, H. Brune, and 
K. Kern, Phys. Rev. B 59, 15926 (1999) 

D. Jacob, M. Soriano, and J. J. Palacios, Phys. Rev. B 
88, 134417 (2013) 

°° J. Kiigel, M. Karolak, J. Senkpiel, P.-J. Hsu, G. Sangio- 
vanni, and M. Bode, Nano Letters 14, 3895 (2014) 

S. Karan, D. Jacob, M. Karolak, C. Hamann, Y. Wang, 
A. Weismann, A. 1. Lichtenstein, and R. Berndt, Phys. 
Rev. Lett. 115, 016802 (2015) 




